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Abstract 

Using real-space renormalisation techniques we analyse the Ising model on 
a Sierpihski gasket with anisotropic microscopic couplings, and observe a 
restoration of isotropy on macroscopic scales. In particular, via use of a 
decimation procedure directly on the fractal lattice, we calculate explicitly 
the exponential anisotropy decay coefficients near the isotropic regime for 
both ferromagnetic and antiferromagnetic systems. The results suggest the 
universality of the phenomenon in lattice field theories on fractals. 

I. INTRODUCTION 

Restoration of macroscopic isotropy in a system which has an underlying microscopic 
anisotropy is a phenomenon which, while absent in uniform media, appears almost universal 
in fractals. Recently, Barlow et al. |T[ considered this phenomenon in the form of resistor 
networks. In particular, they considered a resistor network constructed on the Sierpihski 
gasket with locally anisotropic resistance elements. Then, by successive use of star-triangle 
transformations they obtained recursion relations for the resistance elements from one length 
scale to another, and observed a restoration of isotropy at the macroscopic scale. The exten- 
sions to other fractal objects were also considered and the results suggest the universality of 
the phenomenon at least on exactly self-similar fractals. A question which naturally arises 
in this regard is whether the same mechanism is present for systems in statistical mechanics, 
or quantum field theories. 
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The statistical mechanics of spin systems on fractal lattices, on the other hand, has been 
a subject of considerable interest for which a number of models has been studied. In par- 
ticular, for exactly self-similar fractals with a finite degree of ramification, exact real-space 
renormalisation techniques can be pursued for simple spin models, and as a consequence the 
Ising spin model with isotropic couplings has been solved exactly on various fractals. See, for 
example, Gefen et. al. HH, Luscombe and Desai [§] or Stosic et. al. || for further details. 
The physical motivation behind these investigations is in part that the results obtained on 
fractals may have some relevance for real random systems, or crystals with defects. Also, 
it is interesting to note that this analysis indicates properties of spin systems in noninte- 
gral dimensions which are not apparent in more formal analyses such as the e-expansion, 
and also provides the possibility for considering how this physical behaviour crosses over to 
uniform media ||. 

In the present Letter, we study the Ising model on a Sierpihski gasket with locally 
anisotropic coupling configurations. As a comparison to the formulation presented in [|T|, 
we first analyse this model with a successive use of star-triangle transformations and one- 
dimensional decimation, in order to study the level of anisotropy at larger length scales. 
The results obtained show enormously rapid recovery of isotropy at larger scales. However, 
unlike the case of resistance networks, the use of the star-triangle transformation in the 
present context is not necessarily advantageous in the sense that frustrated spin systems 
cannot be studied within this methodology. Nevertheless, as we noted above, in the case of 
the Ising model, exact real-space renormalisation transformations can also be performed in 
the frustrated system (see for example Stinchcombe || and Grillon & Brady Moreira @), 
and as a consequence we are able to study the large scale structure for both ferromagnetic 
and antiferromagnetic Ising models. In both cases, we observe recovery of isotropy at an 
exponential rate. More specifically, the technique used allows us to rigorously prove the 
convergence of the system to isotropy and to investigate the rate at which homogeneity is 
restored. In particular in the near isotropic regime the scaling to isotropy is exponential 
with a coefficient given by p = 1 ± tanh2 with the upper(lower) sign corresponding to the 
f erro ( ant if erro) - magnet ic system . 



II. RESULTS FROM STAR-TRIANGLE TRANSFORMATIONS 

The model is constructed as follows. We embed the Sierpihski gasket in two Euclidean 
dimensions and place Ising spins o,i = ±1 at each vertex. Geometrically the Ising model is 
then effectively carried by a space of Hausdorff dimension d = In 3/ In 2 = 1.585 with a finite 
degree of ramification R ma x = 4. We consider the standard nearest neighbour Hamiltonian 
given by 

H = -J2 K ijWj , (!) 

(«) 

where = Jij/kT are anisotropic ferromagnetic couplings which depend upon the orien- 
tation of the spins cij and <7j, and the summation is taken over all the nearest neighbour 
sites on the fractal lattice. 
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FIGURES 




FIG. 1. A stage in the renormalisation group procedure for the Ising model on the SG is 
indicated. The interior spins ^1,^2,^3 ar e integrated out resulting in course grained couplings 
L n +i and K n +i. This is achieved in (A): via a combination of star-triangle transformations and 
ID decimation, the dashed lines indicating the transformation used at each stage (the various 
couplings being defined in the text); while in (B): direct decimation on the fractal lattice is used. 



Motivated by the technique used in ]I| we can construct a real-space renormalisation group 
procedure as indicated in Fig. 1 [part (A)] where after the n th iteration of the process we 
denote the horizontal coupling by K^° r = K n , and the two 'vertical' couplings as K" ertl = 



K 



vert2 
n 



r n K n so that r n represents the degree of anisotropy after the n iteration 



The process involves the use of an alternate sequence of star-triangle transformations 
and one-dimensional decimation so that if we define the following functions 



F 1 (a,^) 
F 2 (a,P) 
F 3 (a,f3) 
F A {a,f3) 



|arccosh [e 2a cosh 2/3 



arctanh 



I In 
±ln 



e 2a sinh 2/3 



sinh(2Fi(o, / 3)) 



a±f3 
a-/3\ 
' 2a+f3 ~ 
2a-f3 

cosh 2 2(5 



(2) 



sinh 2 2/3 tanh 2 a 



then the couplings at each stage of the process indicated in Fig. 1 are given by 
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' g = Fi(K n , L n ), h = F 2 (K n , L n ) 

M = F 3 (g,h), N = F 3 (h,h) 

< W = Fi(iV,M), V = F 2 (N,M) (3) 

x = F 3 (g,V), y = F 3 (h,W) 

K n+1 = F 4 (x,y), L n+1 = F 5 (x,y). 

This procedure may be iterated numerically, and with the initial conditions K = 0.01, 
L = 1, i.e. an initial microscopic anisotropy level of r = 100, the results we obtain are 
shown in Fig. 2. The simulation indicates a rapid restoration of isotropy although the actual 
structure of the anisotropy oscillates after each iteration, due to the fact that an odd number 
of star-triangle transformations are used during each iteration process, as indicated in Fig. 1. 

One can easily observe, from the form of the functions F 1 and F 2 , that a 'cut' along 
the negative real axis of the coupling space prevents the consideration of negative couplings 
and therefore a frustrated system cannot be treated via this procedure. In other words, if 
the initial couplings for each triangle are all negative (antiferromagnetic), then the resulting 
couplings on the 'star' become pure imaginary. 




n 



FIG. 2. The level of anisotropy r n plotted against the order of the coarse-graining iteration of 
type (A). Starting from initial ratio vq = 100, we recover r ~ 0(1) after only three iterations. 

Physically, this reflects the fact that, while on the triangular lattice there are frustrations 
associated with each plaquette, these frustrations cannot be realised for any real-valued 
couplings on the 'star', or equivalently, a hexagonal lattice. 

III. FERROMAGNETIC CASE 

In order to allow us to analyse more general situations, and also consider the possibil- 
ity of obtaining analytic results for the asymptotic behaviour, we now consider a different 
method, based on the direct use of real-space renormalisation transformations. As men- 
tioned above, due to the finite degree of ramification, exact renormalisation group analysis 
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is easily accomplished analytically via decimation, as reported in |2|-Q. A course graining 
transformation of this form is represented in part (B) of Fig. 1. 

Consideration of the partition function allows the recursion relations for the couplings L 
and K to be determined as follows, 



, /v 2 cosh 2K n + cosh 2L n + e 4Kn cosh 6L n 



2cosh 2 2L n (l + e 4K « cosh 2L n ) V ; 



4i 2 cosh 2i^ n + cosh 2L n + e 4Kn cosh 6L n 

, (5) 



2 cosh 2K n (l + e 4i ^ cosh 2L n ) 
which, in the isotropic limit L n = K n , clearly reduce to the result of Gefen et al. 



e 4K " +1 = - ^ ■ o ' - ( 6 ) 



e 8i<"„ _ e 4X n _|_ 4 
3 4X„ _|_ 3 



After n iterations of the renormalisation group transformation the level of anisotropy in the 
interaction is again characterised by r n = L n /K n . Since we are concerned only with this 
ratio of the couplings we can rescale both the horizontal and 'vertical' couplings after each 
iteration by an arbitrary factor. A convenient choice of this factor proves to be 1/K n , and 
as a consequence the couplings at each order now take on the values K n = 1 and L n — r n . 
Clearly this rescaling has no effect on the observable behaviour of the system with regard to 
changes in the level of anisotropy. However, the benefit we gain by using this technique is 
that by taking the ratio of (Q) and (Bp we can construct a recursion relation purely for the 
level of anisotropy in the system given by 

Before we analyse this relation we should first consider the possible initial conditions, i.e., 
the initial choices one can make for r . Ignoring the trivial isotropic case there are clearly 
two; either r > 1 or r < 1. With reference to Fig.l we see that these conditions are not 
equivalent. The first corresponds physically to the situation where couplings in two of the 
directions on the lattice are larger than the third, while the latter case corresponds to the 
situation where two of the couplings are weaker than the third. To avoid confusion we shall 
always differentiate these particular initial conditions and refer to each by the corresponding 
values of r . 

The structure of the relation ([?]) indicates that it possesses a two element fixed point set 
{1, oo}, with the first point r* = 1, corresponding to isotropy. In order to determine which 
of these fixed points is stable we analyse the recursion relation for each of the possible initial 
conditions. Making use of the trivial inequality, (a + b)/(c + d) < a/c for any a, b,c,d G R + , 
we find 

1 e 2r " 

r n+ i < 1 + - ln — — = r n , for r n > 1, (8) 



2"' e 2 

1, e 2r " 
- In — - 

2 e 2 



r n+ i > 1 + - ln —j- = r n , for r n < 1. (9) 
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FIG. 3. The deviation from isotropy R n of the interaction plotted against the the renormalisa- 
tion group iteration n in the ferromagnetic case. The upper plot uses the initial condition ro = 10 
while the lower curve uses ro = 0.1. 

Thus the deviation from isotropy at each order 5 n = \r n — 1| satisfies the inequality 5 n+ i 
unless r n = 1 and thus r* = 1 is the only stable fixed point of the recursion relation and, as 
a consequence, restoration of isotropy on macroscopic scales is assured. 

The actual behaviour of the system as one carries out the course-graining transformations 
may be observed by iterating (|?]) numerically. Specifically, in Fig. 3 we present the scaled 
deviation from isotropy 

(10) 

for examples of the possible initial microscopic conditions r = 10 and r = 0.1. It is clear 
that the presence of the unstable fixed point at r* = oo will only affect the restoration of 
isotropy in the first of these cases. With the initial condition ro = 0.1, recovery of isotropy 
is smooth and quite rapid, although approximately an order of magnitude slower than the 
corresponding result using star-triangle transformations (see Section 2). For later reference 
it is useful to characterise this regime as having 'positive curvature' (i.e., r"(n) > where we 
symbolically treat r n as a smooth function). In contrast to the case r = 0.1, with the initial 
condition ro = 10 restoration of isotropy is quite slow for several iterations and appears to 
be significantly affected by the fixed point at r* = oo. We characterise this regime as having 
'negative curvature' (symbolically r"(n) < 0). We observe from Fig. 3 that after sufficient 
iterations there is a crossover from this negative curvature regime to the positive curvature 
regime, corresponding to dominance of the r* = 1 fixed point, at approximately r n « 1.3. 
As a consequence we can conclude that the negative curvature regime is only present if the 
initial configuration is such that ro > 1.3. 

As we have mentioned above, the phenomenon of macroscopic restoration of isotropy 
appears to be an almost universal phenomenon in fractal media. In order to test this 
conjecture in a more quantitative manner it would be helpful to obtain a measure of the 




R„ 



1 

1 

r„ 
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scaling near the isotropic fixed point r* = 1. A numerical study of the recursion relation 
(0) in this regime, for an arbitrary choice of initial conditions, indicates that the rate of 
restoration of isotropy appears to be approximately exponential. That this is indeed the 
case may be verified analytically as follows. If we treat r n as a function, r = r(n), we 
have in general r(n + dn) = r(n) + r'(n)dn where the derivative is taken with respect to 
n. In the regime where r n ~ 1 the variation in r between iterations is very small and 
thus we may approximate the recursion relation (|?p by the truncated Taylor series, i.e., 
r(n + 1) m r(n) + r'(n). As a consequence, we obtain the differential equation 

dr cosh2r 

d^ + r -" 1 ^hJ- 1 = - 

which, for arbitrary initial conditions, is valid in the regime near isotropy where r is only 
weakly dependent on n. Recalling the definition of 8 = \r — 1|, in the limit 5 — > where 
the differential equation is valid, we may expand the nonlinear term in ([11]) to obtain the 
following simple differential equation for S, 

i r 

— + (1 + tanh2)5 = 0, (12) 
dn 

with the solution given by 

6(n) = exp(— pn), (13) 

where the decay coefficient for the level of anisotropy near the isotropic regime is given by 
p = 1 + tanh2 « 1.96. 



IV. ANTIFERROMAGNETIC CASE 

A significant advantage of the use of direct decimation on the fractal lattice is that, in 
contrast to the results obtained using star-triangle transformations, there are no 'cuts' in 
the recursion relations (f|) and @ for negative values of the couplings. This allows us to 
adapt the technique to consider a frustrated system with negative couplings K and L . 

By rescaling the couplings after each renormalisation group iteration, as discussed in 
Section 3, we can derive a recursion relation for the level of anisotropy in a manner analogous 
to that described for the ferromagnetic case, obtaining 

1 , / cosh2r n \ A . 

r "« = 1 -2" 1 l^i^j- (14) 

The possible initial conditions for this relation again fall into the same two classes, r > 1 
or r < 1, which were relevant in the previous discussion. However, for this model the fixed 
point set for the recursion relation (113) reduces to a single point {1} and from this fact 
alone one might expect that the two a priori different initial conditions will not lead to 
significantly different behaviour of the system under renormalisation. In other words, the 
existence of only one fixed point should lead to a single form of scaling behaviour. 

Before we test this conjecture numerically we can verify the stability of the isotropic fixed 
point via arguments similar to those presented in Section 3. Treating the possible initial 
conditions separately we obtain 
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r n+1 > 1 - - In 



1 , e 2rn 



r n+ i < 1 - - In 



1 e 2r„ 



2-r„ 



2 /v. 



2 - g2 - -n, 



for r n > 1, 
for r„ < 1. 



(15) 
(16) 



Thus the behaviour of r n in the frustrated case is oscillatory about r = 1 between iterations, 
which clearly accounts for the loss of the fixed point at oo in this case. From these relations 
one may again derive the inequality 5 n+ \ < 5 n unless r n = 1 and thus r* = 1 is again a stable 
fixed point of the recursion relation and the system will tend to isotropy on macroscopic 
scales. It is of interest to note that this is similar to the behaviour observed using star- 
triangle transformations in the ferromagnetic system. 

As anticipated, numerical iteration of the recursion relation (|TJj) indicates that for the 
frustrated case there is no qualitative difference in the behaviour of the deviation R n from 
isotropy due to a different choice of initial conditions. Specifically, in Fig. 4 the scaled devia- 
tion R n is plotted for the two initial configurations ro = 0.1 and ro = 10 indicating explicitly 
that the behaviour is qualitatively the same in each case. This is again a consequence of the 
loss of the second fixed point at r = oo. 
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FIG. 4. The deviation from isotropy R n of the frustrated interaction plotted against the the 
renormalisation group iteration n. The upper plot uses the initial condition ro = 10 while the lower 
curve uses ro = 0.1. 

Note that in this case, due to the fact that the oscillatory behaviour is not quite symmetric 
about r = 1, we plot every second iteration for clarity. 

The physical consequence of the loss of the second fixed point at infinity can easily be 
observed from Figures 3 and 4. In the situation where at microscopic scales we have ro >> 1, 
or in physical terms where the interaction of the Ising spins in one direction on the lattice is 
much weaker than than the other two, the restoration of isotropy in the frustrated system 
is significantly more rapid than the standard ferromagnetic system. Thus in this case the 
frustration appears to relieve the anisotropy on macroscopic scales. 

However, if we consider the regime near isotropy we obtain a rather different conclusion. 
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In particular if we analyse the scaling of the 'function' r(n) in a manner analogous to that 
presented in Section 3 we obtain a similar differential equation to (|TT1), that is, 



dr , cosh2r . . 

^ + r + hl ^hT- 1 = °- < 17 » 

The deviation 5j = \r — 1| in this case also satisfies the scaling relation (|T3"|), with a smaller 
decay coefficient p — 1 — tanh2 « 0.036. Thus it appears that in this regime the frustration 
inhibits homogenisation. As a consequence, one can envisage that at least in the situation 
described by an initial condition r >> 1 the presence of a frustration has a markedly 
different effect on the level of isotropy at different length scales when compared to the 
ferromagnetic system. 



V. DISCUSSION 

In conclusion, we have examined the phenomenon of restoration of macroscopic isotropy 
for a simple interacting system, a microscopically anisotropic Ising model on a finitely ram- 
ified fractal, the Sierpihski gasket. This is an interesting example of a phenomenon that is 
absent in uniform media. 

The results for the Ising spin model shows an exponential recovery of isotropy at large 
scales, for both ordered (ferromagnetic) and disordered (antiferromagnetic) systems. In par- 
ticular, it is interesting to note that for disordered systems, the phenomenon is independent 
of the initial conditions, while for ordered systems the behaviour depends sensitively on 
the initial conditions. Furthermore, we have also observed that, in the strongly frustrated 
regime, the restoration of isotropy is considerably slower than the corresponding unfrustrated 
cases. 

The above results are also suggestive that further studies of models on spaces which are 
explicit geometric objects of nonintegral Hausdorff dimensions could uncover other subtle 
characteristics which may be of relevance to various physical systems. It could also be 
conjectured that the scaling behaviour which was analytically determined near isotropy 
may in fact be universal in such magnetic systems. 
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sions. The financial support of D.B. by PPARC and A.R by the Commonwealth Scholarship 
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